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Abstract 



The Luttinger liquid in which the concentration of electrons varies ran- 



J> ■ domly with coordinate is considered. We study the fluctuations of the tunnel 

m _ 

conductance, caused by the randomness in the concentration. If the concen- 
tration changes slowly on the scale of the Fermi wavelength, its prime role 
ON ■ reduces to the scattering of the plasmon waves, propagating along the system. 



As a result of such a scattering, plasmons get localized. We show that the lo- 



e 

*^ • calization length, l u , of a plasmon with frequency u is inverse proportional to 

C3 . 

the square of the interaction strength and changes with frequency as l w oc to . 



If the relative variation of the concentration is small, the randomness-induced 

correction to the tunnel conductance, 8G(V), where V is the applied bias, can 

be expressed through the spectral characteristics of the localized plasmons. 

The magnitude of the correction, (5G 2 ) 1 / 2 /G, increases with V as y/V . The 

typical period of the fluctuations in 5G(V) is of the order of V. At a fixed V, 

the correlator of 5G at different points of the liquid falls off with distance as a 

power law and oscillates with the period which is one half of the wavelength 

of a plasmon with frequency u = eV/fi. 
PACS Numbers: 73.20.Dx, 73.40.Gk 
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I. INTRODUCTION 



It is well-known that the density of states in a pure one-dimensional interacting electron 
gas (Luttinger liquid) vanishes as a power law in the vicinity of the Fermi levefi 
where the exponent is determined by the interaction strength. Such a behavior should 
reveal itself in the dependence of the differential tunnel conductance on the applied bias: 
G{V) = dl/dV oc v(eV). It is obvious that the presence of a disorder would perturb the 
local value of the density of states and, thus, cause some random correction, 5G(V), to the 
conductance. Then the relevant questions are: 

i) What is the typical magnitude of SG for a given realization of the disorder? 

ii) How the values of SG at different voltages are correlated? 

iii) How the correlation between the values of 5G at the same voltage, but at different 
points of the liquid, falls off with increasing separation between the points? 

These questions are addressed in the present paper. We will consider the case of a 
smooth disorder. Namely, we will assume that the correlation radius is much larger than 
the Fermi wavelength. This assumption simplifies the problem drastically, since it permits 
one to neglect the backward scattering of electrons^ and, hence, to view the disorder as a 
random variation of the electron concentration with coordinate. The key to understanding 
the role of the disorder is provided by the bosonization procedure!!^, which allows to de- 
scribe the low-energy excitations of the system in terms of bosonic excitations (plasmons), 
propagating along the liquid. Since the velocity of a plasmon depends on the concentration 
of electrons, the spatial variation of the concentration would give rise to the backscattering 
of plasmons. In other words, in the presense of a disorder, a plasmon with frequency uj 
acquires a finite mean free path l^. It is important to note that a smooth disorder, for 
which the backscattering of electrons at the Fermi level is suppressed, might, ultimately, be 
quite efficient in backscattering of plasmons with wavelengths of the order of the correlation 
radius. 

For non-interacting electrons in one dimension it is established that even a weak disorder 
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leads to the localization of all eigenstatesIH. Then the mean free path acquires the meaning of 
the localization radius. The same conclusion applies, certainly, to plasmons. The difference 
is, however, that an electron becomes more and more localized as its energy decreases, 
whereas for plasmons the situation is the opposite: the lower is the frequency, the weaker 
is the localization. Since the dispersion law of a plasmon is linear: uj = sk, where s is the 
sound velocity, the density of the plasmon states is frequency-independent. The reason for 
suppression of the localization for the low-frequency plasmons is that the matrix element of 
the baskscattering vanishes in the limit of long wavelengths. In particular, we show that for 
small us the frequency dependence of the mean free path is l w oc oj~ 2 and, correspondingly, 
the product kl^, which determines the localization strength, increases as 1/u when u goes 
to zero. We also show that l w turns to infinity when the interactions are switched off, thus 
revealing that the backscattering of plasmons is possible only due to interactions. 

If the relative variation of the concentration is small, the disorder-induced correction 
to v(u) can be expressed in terms of the spectral characteristics of the localized plasmons 
(such as local density of states). On the other hand, these characteristics were the subject 
of detailed studies in application to the localized electrons®!. By utilizing the approach 



developed in Refs. we calculate the variance (SG 2 ) 1 ^ 2 , and the two-point correlator 



5G(x)5G(0). We find that the ratio (5G 2 ) 1 / 2 /G increases with voltage as yV and that the 
characteristic "period" of change of SG(V) with voltage is of the order of V. We also find 
that the two-point correlator falls off with x as a power law and oscillates with the period 
5x = 7csh/eV, which is one half of the wavelength of a plasmon with frequency u = eV/H. 
The paper is organized as follows. In the next section the formula for the mean free 



path of a plasmon is derived. In Section |TTTj we calculate the correlator of fluctuations of the 
tunnel conductance. Section |TV] concludes the paper. 
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II. MEAN FREE PATH OF A PLASMON 



The Hamiltonian of a Luttinger liquid with concentration of electrons, n(x), being a 
function of coordinate has the formSini 



n 



dx 



p 2 (x) 



1/ 7T 2 h 2 , v 

2 V m 



d(nu) 



(1) 



_2mn(x) 2 V m v ' ) V dx 

where u and p(x) are, correspondingly, the displacement and conjugate momentum 
([it(x),p(x')} = ih5(x — x')); Vq is the effective interaction strength: Vq — J dxV(x). The 
Hamiltonian can be reduced to a system of harmonic oscillators 



^ = E 



\V\ mQ 2 . ■ 
2m + 2 ^ 



(2) 



by means of the following transformation 

1 



u[x) 



E 



n{x] 



^y.{x)Q„ , p(x) = E V "(^^(x)^ 



(3) 



Here $ M are the eigenf unctions of the operator D which is defined as 
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D<5>, 



n(x 



A. 

dx 



m \ m I J dx 



n(x)§. 



(4) 



The eigenvalues of D determine the frequencies, Q^, of the oscillators. If the concentration 
is constant (n(x) = no), the solutions of ([|) are the plane waves 

with a linear spectrum f2 M = 0° = s/c^ (L is the normalization length). Substituting ([5]) 
into (H), we get the standard expression for the sound velocity^ 



s = V F \ 1 



V a 



(6) 



where Vp — irhn /m is the Fermi velocity. 

Assuming that the relative variation of the concentration is small, \n(x) ~uq\ <C uq, the 
expression for the mean free time, t w , for a plasmon with frequency u is given by the golden 
rule 



1 



71 



UJ 



E|(^ 0) 



D 



$i 0) )^° )2 -^ 2 ) 



(7) 



where is the plane wave with = lo/s and (• • •) stands for the averaging over the 
fluctuations of n(x). The factor 1/u in (0) appears since we define the mean free time as 
1/t w = ImQ^ = Imf2 2 /2cG>. Keeping only the first order terms in the difference n(x) — no, 
we get the following expression for the matrix element 



dx 2 



n 



G?X — TT-o) 



2 " dx 2 2 " dx 2 dx rfx 



The energy conservation, insured by the 5-function in (|7|), requires that k v = —uj/s 
(backscattering) . Then we have $| 1 - ) * = $[ / °- ) and, as it can be easily seen, the second 
term in the matrix element vanishes identically. As a result, the matrix element takes the 
form 



<*(o> £) *?>) 



mL \ s 



dx (n — n )e 



2iux/s 



(9) 



We see that the matrix element for backscattering is proportional to the interaction strength, 
which reflects the fact that this process is possible only due to interactions. The mean free 
path, l u , defined as l w = sr^, can be found after substituting (§) into (|7|) 

1_ f s 2 - v] 
L ~ V s 2 



(5n) 2 /uj\ 2 / + x 



2n 2 



dxK ( — 

oo V R r 



(10) 



where we have introduced the correlator of the fluctuations of nix) 



(n(x) — no)(n(x') — n ) = (5n) 2 K 



Rr. 



(11) 



Here R c is the correlation radius and (5n) 2 is the mean square fluctuation of n(x), so that 
K(0) = I. 

It is instructive to rewrite fllOD for the product kl u = ul^/s, which measures the effective 
scattering strength. Assuming the gaussian form of the correlator, K[z) = e~ z2 , we obtain 



(kl u ) 



-i 



7T 



(dn) 2 
~2n[ 



loR, 



c e -u 2 B£/s 3 



(12) 



Note that the last factor in ([12]) is a function of the argument uR c / s and this function does 
not exceed unity. Thus, the product kl w is large for any u. Obviously, the backscattering is 
ineffective when ujR c /s > 1, i.e. when the correlation radius R c exceeds the wave length of 
the plasmon. In the opposite limit, ujR c /s < 1, the mean free path behaves as ~ oj~ 2 . As 
it was mentioned in the Introduction, the real meaning of l w is the localization length of a 
plasmon with frequency to. The fact that kl^ > 1 allows one to apply to localized plasmons, 
the description developed earlier for localized electrons^ just in this limit. This is done in 
the next section. 



III. CORRELATOR OF THE FLUCTUATIONS OF THE TUNNEL 

CONDUCTANCE 

Once the Hamiltonian is diagonalized, the derivation of the formula for the density of 
states becomes standardly. The operator ^(xo), creating an electron at point xq, can be 
presented in the formlllfl 

^+(x ) = exp E t>(zo)^ + PrMQp) j , (13) 
where the coefficients a M and (3^ are defined as 



a li (x )= dx J!— , Ph(xq) = 7v^n(x ) $ M (x ). (14) 



' x ° yn(x) 

The calculation of the Green function (\l/(£)\l/t(0)) does not differ from that for a pure 
Luttinger liquid and leads to the following expression for the density of states 

v{uj) = -Re / dte iu}t (V(tW{0)) = -Re / dte iut e~ w{t \ (15) 

7T JO 7T JO 

where W is the sum over eigenmodes 



«=?T 1+ S:I( 1 -^ (16) 



If n(x) is constant so that the eigenfunctions are given by (B|), Eq.(^5j) reproduces the known 
result@ for the density of states: vq(u) oc uj k , with 



1 / S Vp 



1 + 



Vo 



K = — \ 1 — I - 1 = ; 27ThV ^ ~ 1 

2\v F s J Jl + -g>- 



(17) 



irhvp 

Indeed, for plane waves we have \a^\ 2 = l/k 2 Ln , |/3 M | 2 = n 2 n /L and we get the following 
expression for the function W 



W (t) = l(^ + V -f) J™ f (1 - e™)e-™ = ~ ( - + - ] In 



2 \vf s J \ r 



;is) 



where r is the cutoff parameter. Substituting this expression into (|15|) leads to Ql7|). 

To find the disorder-induced correction to vq(uj), we treat the difference W(t) — Wo(t) as 
a perturbation and expand the exponent in (O) to the first order. This gives 



Su 



-Re/ dfe*" 4 ' ° 



(19) 



7T Jo vro + zst, 

To calculate the integral over t, it is convenient to present the denominator in (|T9"| ) in the 
form 

1 1 



dzz K e~ z( ™\ 



(20) 



(r + ist) K+1 T(k+1)Jo 
After substituting ( PU[ ) into (|T9"D, both integrations, over z and i can be easily carried out, 
and one obtains 



Su(xq, uj) 



1 



r 



r(« + i) v s 



K+l 



E 



m\a^\ 2 ^^ ^1/3^1 



2ft 



2mfi, 



(a; - fi. 



r /s ^ (a; _ sg )4 (2i) 

JO Q J 

by integrating 



Note that the constant a u can be expressed through the derivative ^f-f 
the equation ([|) from s to oo. Then, one gets 

s 2 1 /ci$„ N 



01 r 



Q 2 ^/hq \ dx 



(22) 



To study the correlation properties of 5u, we introduce the following local densities 



pi(e,x ) =^|$ M (x )| 2 5(e-fi M ), 



p 2 (e,x ) = s 2 J2tv 



rf$ M (x ) 



5(e - n, 



(23) 
(24) 



The average values of these densities are equal 



Pi = P2 



1 



7XS 



(25) 



Then the correction 8v can be rewritten in terms of the fluctuations 5p± = p\ — 1/tts and 
Sp2 = P2 — 1/tts 

1 









i T Jo 



-{u-u x ) K 



—Spifaxo) + —5p 2 {uj 1 ,Xo) 
s v F 



• (26) 



Correspondingly, the correlator of 5v at different points is expressed through the correlators 
of the fluctuations Spi,Sp 2 

2 



8u(x\, co)5p(x 2 , uj) 



1 



T(k + 1) V s ) 2 
u duj\ du>2 



X 



10 IQ\ JO U) 2 

where the correlators Fi, F 2 and F\ 2 are defined as 



(u - u^u - oo 2 r -f Fi + — F 2 + 2F 12 , 



(27) 



Ft = 5p 1 (x 1 ,uj 1 )5p 1 (x 2 ,uj 2 ) , 
F 2 = Sp 2 (xi,uj 1 )Sp 2 (x 2 ,uj 2 ) , 
F 12 = 5p 1 (x 1 ,uj 1 )6p 2 (x 2 ,uj 2 ) . 



(28) 
(29) 
(30) 



The correlator F\ is, in fact, the correlator of the fluctuations of the local density. In 
application to localized electrons, it was studied in Refs. §,[| using, correspondingly, the 
Berezinskii techniquei and the technique developed by Berezinskii and Gor'kovi. In the 
limit cuit Wi 1, co> 2 t W2 1, the correlator is nonzero only if the difference U\ — uj 2 is small 
enough: \uj\ — uj 2 \ ~ l/r Wl . For distances z = \x\ — x 2 \ <C Li,L 2) the expression for F\, 
obtained in Ref. |8|, reads 



3vrV 



2r, 



T 6(e) (3 - 2 sin 2 (^)) + 2 sin 2 (^) (C(s) - 1) 



an 



(31) 



where e = io\ — uj 2 is the difference between the two frequencies and the function C(u) is 
given by 
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C(s) = r dy ™^T y) = roo rfge - g Q (32) 

h y + 1 jo q z + 

As it was noted in Refs. |9|,j8|,[13|, at z = and e 7^ 0, we have -F\ = 0, i.e. the correlation is 
absent. 

It is easy to establish that the correlator F 2 is equal to F\. Concerning the correlator F± 2 , 
we did not find the expression for this correlator in the literature. So, we have calculated it 
using the Berezinskii-Gor'kov technique and obtained the following expression 

1 
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71 5(e)(l + 2sin 2 (^))+2cos 2 (^)(C(,) 



., , , o \ ' - ■ x , / - i / ' < ^ ( 33 ) 

It is seen that in contrast to F\ , the correlator F\ 2 is non-zero at z = and e 7^ For finite 
e, both correlators are proportional to (C(e) — 1) and decay with increasing e as (st Ui )~ 2 . 

To calculate the double integral in (0), we make use of the fact that the correlators 
Fi,F 2 and F 12 are sharp functions of e, i.e. the major contribution to the integral comes 
from the domain | — 0^2 1 ~ l/ r wi ^ uj 1 ,uj 2 . This allows to put uj 2 = u)\ in all other factors 
and to extend the integration over £ = io\ — uj 2 to (—00,00). Noting that the integral of 
(1 — C(e)) is equal to k/t u1 , we obtain 

, / v% s 2 „ „ \ 1 (a v F \ 2 f2uj 1 z\ , . 

de[^F 1 + -^F 2 + 2F 12 ) = —g— (_ _ -Z) cos(-^). (34) 



00 



Note that the coefficients in front of Fi , _F 2 and F 12 in ( |2"7| ) have combined into the factor 
(s/vf — vf/s) 2 which is proportional to V 2 at small Vq. Using (|3"4j), the expression for the 



correlator 5v(xi,uj)5is(x 2 ,uj) takes the form 

— — r vr / s v F \ 2 [ 1 /r \ K+1 l 2 dwi , , 2K /2cui2\ 

5z/ x 2 ,w = — — — — / — — (w-ui) cos 

8 \vf sj T{k + 1)\sJ Jo u{t w , \ s J 



(35) 

Apparently, the factor 1/uj 2 in ( p5|) diverges at small U\. However, this divergency is com- 
pensated by the frequency dependence of the mean free time, which at small U\ behaves as 
1/u; 2 . Substituting r Wl in (^) and introducing the new variable w = ui/w, we get the final 
result 
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5v(xx, u)5u(x 2 , to) 



71 



3/2 



16(2k + 1) 



r(« + 1) V s 



r 



K+l 



(Sn) 2 sR c u 2K+1 



n 2 Q 



R(z), 



(36) 



where the function R is defined as 



R(z) = (2k + 1) f 1 dw(l - wf R e 
Jo 



2k -wH^) 2 



COS (WZ) 



(37) 



Here, z is the dimensionless distance 



Z = [X2 - Xl) ■ 

s 



(38) 



The function R is defined in such a way that, in the only interesting limit ujR c /s 1, it 
turns to 1 at z — 0. This function determines the coordinate dependence of the correlator 
( p7[) and, consequently, the coordinate dependence of the correlator of the fluctuations of 
the tunnel conductance 



5G(x)5G(0) _ /2eV ^ 



(39) 



(<5G(0)) 2 " V sh 

The function R is ploted in Fig.l . With increasing distance, it falls off and oscillates. The 
asymptotic behavior of R(z) at z 3> 1 is as follows 



I}/ ~\ T(2k + 2) . 
R\z) _„ - sm(.: - ttk 1 



g2n+l 



(40) 



The spatial period of the oscillations of R is Sx = Txsh/eV and it decreases with increasing 
voltage. It is to be emphasized that such a behavior of the correlator of 5G (slow decay and 
oscillations) is entirely due to interactions. In the absence of the interactions, the correlator 
( |55D would simply reproduce the correlator ( |Tl"D of the fluctuations of n(x), i.e. it would 
decay monotonously at distance ~ R c which is much smaller than Sx. 

At z — 0, the formula fl2"T|) defines the variance of Sv. It is convenient to rewrite it for 
the relative magnitude of the fluctuations of SG 

^Jl(^) 1/2 (41) 



((SG) 



2U/2 



7T 



3/4 



G 



4(2re + l) 1 /2 v F 



n 
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Naturally, the variance (5G 2 ) is proportional to the magnitude of the fluctuations of the 
density. Less trivial is that ( f41|) contains two additional factors: (1 — vp/s 2 ) 2 which is 
proportional to Vq at small Vo, and a small factor {ujR c /s) 1 I 2 (note that ([H]) is written in 
the limit uj <C s/R c ). This factor determines the voltage dependence: {{5G) 2 ) 1 ^ 2 / G oc y/V. 
Qualitatively, this dependence can be interpreted as follows. Formula (|26| ) shows that all 
plasmon modes with frequencies smaller than u contribute to 8v(xo,u). Obviously, the 
major contribution comes from plasmons for which the center of localization lies within the 
localization radius from the point xq. Let us consider some frequency strip (ui — A, uj\ + A) 
centered at u)\ < uj and with the width A<a;i. Then, the localization radius for all modes 
within this strip is approximately l m . For a spatial interval of length l u , the typical frequency 
spacing between modes is l/pl Wl ~ 1/t Wi . Correspondingly, the average number of modes 
within the strip is At LUi and, hence, the relative fluctuation of this number is ~ (At^) -1 / 2 . 
The final estimate for the fluctuation of 8v emergies if we set A ~ uj\ ~ uj. Then, one gets 
bv/v = 5G/G ~ {ujt^ 1 / 2 w {kl^ 1 / 2 . Substituting (JT2I) for kl w , we reproduce the voltage 
dependence in fl41[) . However, the estimate and the result of the calculation still differ by 
a factor (1 — vp/s 2 ). This factor originates from the specific details of the strucure of the 
eigenfunctions (see the comment after Eq. fl34j) ), and we cannot interpret it qualitatively. 
Note in conclusion of the section, that since all frequencies from to w give, roughly, 
comparable contributions to 8i/(xq,u), it will change significantly only when the frequency 
doubles. This means that at fixed xo the characteristic period of fluctuations in 5G as a 
function of voltage is of the order of V. 



IV. CONCLUSIONS 

The main result of the present paper is that the randomness in the concentration of 
electrons in the Luttinger liquid causes a random component in the tunnel conductance 
which changes semiperiodically along the liquid, with the period depending on the applied 
bias. This behavior results from the fact that at a fixed bias, V, the frequencies of plas- 
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mons, responsible for the correction to the density of states, are strictly limited by eV/H. 
Correspondingly, the oscillating behavior of the correlator of the fluctuations reflects the 
distribution of density in a plasmon with maximal frequency. If the conductance is studied 
as a function of bias, the disorder would cause fluctuations with characteristic period of the 
order of V. The fluctuations disappear as the wavelength of the plasmon with frequency 
eV/h becomes smaller than the spatial scale of the change of the concentration, R c . There 
is also a limitation from low voltages, imposed by the finite length of the liquid, L. Namely, 
our calculation applies when the mean free path of a plasmon, 1^, is smaller than L. At 
low enough frequencies l u exceeds L and the disorder does not play any role. For such 
frequencies the oscillations of the tunnel conductance with voltage have their origin in the 
size quantization of plasmons, the period of oscillations being irhs/eL. These oscillations 
were studied in Ref. |I0| . 

Note finally, that the approach developed in the present paper can be extended to the 
case of the multichannel Luttinger liquidS In the latter case, the scattering of plasmons 
between the channels, caused by a disorder, should be taken into account. 
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FIGURES 

FIG. 1. The correlation function of tunnel conductance is plotted as a function of the dimen- 
sionless distance z for s/vf = 1-2 
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